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Abstract 



We have performed a detailed investigation of geophysical constraints on 
the possible admixture of mirror matter inside the Earth. On the basis of 
the Preliminary Reference Earth Model (PREM) — the "Standard Model" of 
the Earth's interior — we have developed a method which allows one to com- 
pute changes in various quantities characterising the Earth (mass, moment of 
inertia, normal mode frequencies etc.)due to the presence of mirror matter. 
As a result we have been able to obtain for the first time the direct upper 
bounds on the possible concentration of the mirror matter in the Earth. In 
terms of the ratio of the mirror mass to the Earth mass a conservative upper 
bound is 3.8 x 10~ 3 . We then analysed possible mechanisms (such as lunar 
and solar tidal forces, meteorite impacts and earthquakes) of exciting mirror 
matter oscillations around the Earth centre. Such oscillations could manifest 
themselves through global variations of the gravitational acceleration at the 
Earth's surface. We conclude that such variations are too small to be ob- 
served. Our results are valid for other types of hypothetical matter coupled 
to ordinary matter by gravitation only (e.g. the shadow matter of superstring 
theories). 

95.35.+d, 91.35.-x, ll.30.Er, 12.60.-i, 14.80.-j 
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I. INTRODUCTION 



Recently there has been considerable interest in the study of particle physics implications 
of the SuperKamiokande data on atmospheric neutrinos [0] as well as the results of other 
neutrino experiments @ . The main purpose of such investigations has been to understand 
what theories could be responsible for the observed experimental features which give strong 
evidence for large angle neutrino oscillations. One of these theories is the Exact Parity Model 
(EPM) which introduces parity or "mirror" partners for all ordinary particles (except the 
graviton) and thus restores the parity invariance apparently broken by weak interactions 
. The Exact Parity Model predicts pairwise maximal mixing between ordinary and 
mirror neutrinos and provides a basis for interpretation of atmospheric neutrino and solar 
neutrino data M. 

An important question that arises naturally is whether or not the existence of mirror 
particles can lead to other observable consequences. In this work our objective is to find 
constraints on the possible concentration of mirror particles in the Earth. These constraints 
are a necessary step in the search for mirror world manifestations in terrestrial experiments. 
In particular, the presence of mirror matter in the Earth could lead to the regeneration 
of ordinary neutrinos from mirror neutrinos and consequently to the suppression of the 
day-night effect in the neutrino data. 

Two main approaches to our problem are possible. First, one can trace the fate of the 
mirror particles starting from the early Universe epoch through the structure formation 
periods (galaxies, solar system and finally the Earth). Second, we can use geophysical data 
to get a more direct limit on the concentration of mirror matter in the Earth regardless of 
possible cosmological bounds. 

It has been suggested that considerations based on the structure formation theory dis- 
favour a significant presence of mirror matter in the Earth (see Blinnikov and Khlopov in 
H; Kolb, Seckel and Turner in 0). However, as our knowledge of the structure formation 
is still incomplete, it is important to develop a geophysical approach as an independent, 
complementary tool of analysis exploiting the wide and rich variety of observational data 
accumulated in the Earth sciences. 

This approach will be applicable not only to the specific EPM model, but also to any 
other theory predicting the existence of a new world of particles which couples to the or- 
dinary matter only through gravitational interaction. An example is the shadow matter 
characteristic of superstring theories. 

The plan of the paper is as follows. Section 2 summarises the main aspects of the 
"Standard Earth Model" called "Preliminary Reference Earth Model" (PREM). In Section 
3 we analyse the possible effects of mirror matter within the context of PREM. Section 4 
is devoted to the study of dynamical manifestations of mirror matter and comparison with 
gravimetric data. Finally, we present our conclusions in Section 5. 

II. PRELIMINARY REFERENCE EARTH MODEL (PREM) 

The Preliminary Reference Earth Model (PREM) || is a mechanical model of the average 
internal structure of the Earth based, mainly, on the analysis of seismological data. It gives 
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the radial distributions of mechanical properties (such as density, elastic moduli, pressure, 
gravity and others) in the Earth's interior. 

The set of initial data used for constraining the model includes: 

1. astronomic-geodetic data (radius, mass and moment of inertia of the Earth); 

2. data on free oscillations and long-period surface waves (over 1000 eigenfrequencies are 
known) ; 

3. body waves data ( ~ 10 6 arrival times have been registered). 

Let us summarise the main analytical relations used in the construction of the model 
(for more details see e.g. @,]nj). The velocities of the elastic waves are given by 



v p {r) 



K{r) + fMr) 

v s (r) 



Mr) (1) 



where v p is the velocity of longitudinal waves, v s is the velocity of transverse waves, K is 
the bulk modulus (or incompressibility) and p is the shear modulus. These velocities can 
be found from seismological observations as functions of radius r. Measurements of wave 
velocities supply the ratios K/p and pj p. To obtain p independently, the Adams- Williamson 
equation, 



dp Gp(r) 



/ A-Ka 2 p(a)da, (2) 



dr r 2 (t> 2 (r) — |f 2 (r)) Jo 

must be used. This equation expresses the condition of mechanical equilibrium between the 
gravitational attraction and the pressure due to elastic compression. G is Newton's constant, 
and the combination of the squared sound speeds in the denominator is called the seismic 
parameter $: 

*(r) = v 2 p (r) - \vl{r) = ^M. (3) 

Equation @ is valid for a chemically homogeneous layer with adiabatic temperature gradi- 
ent. 

Further, the connections between the density profile and the profiles of pressure and 
gravity are given by 

dP dp . . 



and 



G 

9( r ) = "o / ^a 2 p(a)da. (5) 



r 2 Jo 



III. STATIC CONSTRAINTS 

In this section, we will compute constraints on stationary mirror matter from PREM. A 
later section will consider dynamic manifestations of mirror matter. 
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A. Pedagogical warm-up exercise 



In order to get a feeling for what would go wrong if a substantial amount of mirror matter 
were present in the Earth, we consider a simple but unrealistic scenario first. Suppose the 
Earth is actually two concentric spheres, one ordinary and one mirror. Suppose also, for 
simplicity, that the mass density ratio pi(r) / po( r ) of mirror to ordinary matter is independent 
of radius. This is clearly unrealistic and we will relax this assumption in the next subsection. 

In this case, alterations due to mirror matter of the PREM equations ([I], |3|, £|, |^) can 
be accomodated either by rescaling Newton's constant, 

G -> G' = G(l + ^i) (6) 
Po 

(throughout the paper quantities with index will refer to ordinary matter while those with 
index 1 to mirror matter) or by keeping G fixed and rescaling all of the other quantities used 



While the first procedure (rescaling G) seems to be the simplest, the second procedure 
reveals the physical modifications required in the presence of mirror matter more clearly and 
allows one to use PREM results directly. 

From the formulas (|l|) and the Adams- Williamson equation (fj) we deduce that the 
effective parameters can be defined as 



in PREM. 



Peff = Po + Pi, 



(7) 




(8) 



Peff 



Po(l + -), 



(9) 



K eff = K (l + ^), 



(10) 




(11) 
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With these definitions the form of the Adams- Williamson equation does not change: 

(15) 



d Peff _ Peff9eff 



dr $ e// 

Also, the usual equation for the bulk modulus holds true in terms of the effective quantities: 

dPeff _ Peff 



dP eff Keff' 



(16) 



This demonstrates self-consistency of the rescaling procedure. 

Now, for illustrative purposes, consider the case of a 50%-50% mixture of ordinary and 
mirror matter. The effective values for the density, incompressibility and pressure at the 
Earth centre can be taken directly from PREM data: 

p eff = 13.1 g/cm 3 , K eff = 14.2 Mb, P eff = 3.64 Mb. (17) 

Correspondingly, the values for the ordinary matter (at the centre) are obtained by a factor 
of 2 rescaling: 

p = 6.55 g/cm 3 , K = 7.1 Mb, P = 1.82 Mb. (18) 

With these values, iron is definitely ruled out as the main component of the core and, 
therefore, the Earth as a whole. However, from independent evidence we know that iron 
is one the most abundant elements in the Earth. In addition, no other element (with 
significant abundance) can have the properties required by Eq. flilf ). For instance, silicon 
(at zero pressure) has a lower density than iron, but its incompessibility is greater than 
that of iron. Therefore the 50%-50% mixture of ordinary and mirror matter in the Earth is 
incompatible with the observational data encoded by the PREM model taking into account 
knowledge of terrestrial chemistry. 



B. Modified Adams- Williamson equation 

Consider now the realistic case where the mirror matter density does not follow the 
ordinary density. Then the density of the ordinary matter (indexed by 0) would obey the 
modified Adams- Williamson equation: 



dpo _ Gp 
dr r 2 (v 2 v0 - §w 2 ) Jo 



4:%a 2 (p (a) + pi(a))da. (19) 



Also, we have to require that the total mass of the Earth and the moment of inertia are 
equal to their observed values: 



f R<S 4vra 2 (p (a) + Pl (a))da = M e , (20) 
j o 

■/ a 4 (p (a) + pi(a))da = I, (21) 
J o 



3M ffi Ri Jo 
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where 



M e = 5.974 x 10 24 kg, I = 0.3308. (22) 

Let us start with the simplest case: assume that the mirror matter forms a ball of uniform 
density p\ with a radius equal to the radius of the inner core (Rinner core = 1221/cm « i? e /5). 
Within PREM the density profile in the inner core is given by the following function: 

Pprem{t) = Pprem(0) - q-^, (23) 

where 

Pprem (0) = 13.09 g/cm 3 , q = 8.84 g/cm 3 . (24) 

Therefore, we look for the solution of the modified Adams- Williamson equation in the fol- 
lowing form 

2 

Po{r) = PPREM{0)-e-{q + 5)^-, (25) 
assuming that the mirror density pi is small compared to the Earth central density in PREM 



model Pprem(®)- Solving the system of equations (19, ^5|) we find: 



e « 0.84/)!, 5 « 0.9pi. (26) 

Thus we see that if we wish to add to the Earth some amount of mirror matter then 
consistency with the equilibrium equation requires that the density of ordinary matter be 
decreased (as compared with PREM density) by approximately the same amount. Also, it 
can be shown that if Eq. (^) holds then the constraints due to the Earth mass and the 
moment of inertia are automatically satisfied as long as p\ < 0.77 g/cm 3 . 

Such a decrease of ordinary matter density can be achieved in one of two ways: 

1) by changing the chemical composition of the core so that the new composition has 
lower density than the standard; 

2) by decreasing the pressure so that the density of the standard core composition is 
lowered. 

Before we consider these two possibilities let us review briefly the subject of the standard 



core composition (see e.g. ||10|| ). Information about the chemical composition of the core is 
obtained by comparing the mechanical characteristics given by the PREM model with the 
properties of various substances under high pressure. In this way it has been established 
that the core characteristics are close but not equal to those of iron. There is sufficient 
evidence to conclude that some lighter element should be added to iron in order to satisfy 
the geophysical constraint on the core composition. Sulfur and oxygen appear to be the 
strongest candidates for that role although other elements (such as carbon, nitrogen etc) 
cannot be ruled out at present. As examples, cores containing 6-12% sulfur or 7-8% oxygen 
(by mass) have been proposed as possible compositions consistent with the PREM model. 

In order to accomodate mirror matter according to the first method above, we should 
increase the admixture of light elements as compared with the standard levels. However, 
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such an increase would typically raise the sound velocity in the core (for more details see 
flQfl). Requiring that the sound velocity is equal within accuracy of about 0.3% to the 
observed value (see ||) we can estimate that 

e< 0.18 g/cm 3 . (27) 

Next, in the second method we have to lower the pressure in order to obtain lower density of 
the (ordinary) matter. Lowering the pressure would lead to a decrease of the sound velocity 
and from the same requirement as above we can conclude that 

e < 0.06 g/cm 3 . (28) 

Then, using Eq. ([26|, |27|, [28|) we can obtain the upper bound on the mirror density in the 
inner core 

pi < ^-max{0.18, 0.06} g/cm 3 = 0.21 g/cm 3 . (29) 

Translated to the upper limit on the ratio of the mirror mass to the total mass of the Earth, 
this becomes: 

^< 27X10- (30) 



C. Constraints from free oscillations of the Earth 

Let us consider now the effect of changing the ordinary density on the Earth eigenfre- 
quencies. Using Eq.(A2~A6) of Ref. || and Eq.(41) of Ref. [IT] it can be shown that the 



relative change of period of the spheroidal oSo mode as a result of changing the inner core 
density by e would be 

ST ^ 0.026e 

~T ~ ~Pprem(0Y ( } 

Requiring that this shift of period be less than the fitting accuracy, 0.05%, we obtain a 
bound on the allowed change of the ordinary matter density in the inner core: 

€ —- - < 1.9 x 10" 2 , e < 0.26 g/cm 3 . (32) 

PPREM\y) 

Using Eq. ( p6|) this bound can be translated into a limit on the mirror matter density in the 
inner core: 

< 2.3 x 10- 2 , Pl < 0.3 g/cm 3 . (33) 

Pprem\}j) 

In terms of the total mass of the mirror matter M\ we can rewrite Eq. (|33"D as 

^< 3.8X10- (34) 
Comparing Eqs. (|30|) and ([34]) we see that the difference between these upper limits is not 



very significant, although Eq. ( |30|) is perhaps a less reliable estimate than Eq. ( pl|) because 
of the incomplete knowledge of mechanical properties of various materials at high pressures 
characteristic for the Earth's centre. For these reasons we interpret Eq. fl3"4|) as our final 
conservative upper bound on the mirror matter mass located inside the inner core of the 
Earth. 
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D. Arbitrary radius of the mirror matter ball 

So far we have considered the case when mirror matter is contained completely inside the 
inner core of the Earth. To justify such an assumption one would need to know the detailed 
macroscopic properties of mirror matter (such as equation of state, chemical composition 
etc.); then using the condition of mirror matter equilibrium one could obtain the relation 
between the mirror matter and its radius. If we do not want to rely on such additional 
information than we have to regard the radius of the mirror ball Ri as a free parameter of 
our model. Of course the resulting constraints will be weaker than they could be otherwise. 

For simplicity we will restrict ourselves to the two characteristic values of R\ in addition 
to the case R\ = Rinner core considered before: first, Ri = R ou ter core ~ 0.55i? ffi , and, second, 

R\ = Rlower mantle 

~ 0.89-R®. The second choice is motivated by the fact that for radii larger 
than Ri ower mantle the Adams- Williamson equation is not valid anymore. 

Also, we will assume that the mirror matter has a uniform density pi while the density of 
the ordinary matter differs from its PREM value by a radius-independent correction Spo for 
r < Ri and coincides with the PREM value for r > Ri. Although both these assumptions are 
clearly unrealistic, they nevertheless give us a self-consistent approximation scheme because 
both pi and Spo are small and therefore the effect of their radial dependence would be a 
second-order correction. 

The precise details here depend on the chemical composition, the equation of state and 
other thermodynamic parameters for mirror matter. For instance, if we are given the equa- 
tion of state for the mirror matter then (using some additional simplifying assumptions, 
such as chemical homogeneity, neglecting temperature etc.) we could find the density profile 
of the mirror matter (in particular, its central density) from the equation of mechanical 
equilibrium. However, our goal in this paper was to obtain limits on the mirror matter that 
would be independent of such particularities. 

In any case, our method allows one to calculate the upper limit on the mirror matter 
mass for an arbitrary distribution of mirror matter p\ (r) . 

As before, our constraints on the mass of the mirror matter will be based on 4 pieces of 
information: 

1) mass of the Earth and its coefficient of inertia; 

2) validity of the modified Adams- Williamson equation; 

3) periods of Earth's free oscillations; 

4) velocities of elastic waves. 

Generally speaking, points 1) and 2) tell us that the correction to the ordinary matter 
density should be approximately equal to the density of the mirror matter 5p — Pi- Next, 
using that equality we can compute the shift of the period for the <So mode and then obtain 
the upper limit on the mirror mass. From what follows it will be evident that we do not need 
to find an exact relation between p\ and Spo as it would not significantly change the final 
constraints. Therefore, the validity of the modified Adams- Williamson equation (MAWE 
for short) can be analysed in a simpler manner than that of Sec. B. 

As a criterion of validity of MAWE we can require that the actual mass of the Earth 
(that is, the ordinary mass plus the mirror mass) inside any radius r < R\ should be equal, 
with the accuracy of 1% ||, to the PREM mass of the Earth (within the same radius): 



S 



w = 



M(r) 
M(r) PRE M 



< 1%, (35) 



where 



M(r) = M(r) PRE M + f 4na 2 ( Pl - 5p )da ~ M{r) PREM + ^vrr 3 ( Pl - 5p ). (36) 

J o 3 

It is convenient to introduce the dimensionless ratio 

/ = (37) 

where /% = 5.5 g/cm 3 is the average density of the Earth. In terms of this quantity the 
condition ( pop can be rewritten as 

W =(jt) 3 ( m(1 (S )*/S0.01, (38) 
V-Rffi/ \M{r)p REM J 

for r < Ri. 

Next, we require that the total mass of the Earth equal the observed value with an 
accuracy of 1.3 x 10~ 4 |J which translates into 



RA 3 



x / < 1.3 x 10" 4 . (39) 



Further, we have to demand that the coefficient of inertia of the Earth equals its PREM 
value I prem = 0.3308 with an accuracy of AI / 1 < 3 x 10~ 4 . In terms of / this condition 
reads 

We observe that this condition does not give us an independent constraint because it is 
satisfied automatically as long as inequality ([39]) is fulfilled. 

Finally, we compute the shift of the period for the o5o normal mode of the Earth due to 
non-zero Spo using the same method as in Section C. In the case R\ = Router core we find: 

% « -0.136 (&. - f). (41) 
T VP© / 

From inequality (^) we conclude that 

/<8x 10" 4 . (42) 



With these values of / the criterion (j38"D of MAWE validity is clearly fulfilled for all radii 
r < Ri. Next, requiring that the shift of period be less than the fitting accuracy (0.05%), 
we obtain the upper limit on the mirror matter density 

p x < 0.025 g/cm 3 . (43) 



9 



This translates into the following bound on the mirror matter mass: 



^<7.4xl0" 4 (far R x = 0.55^)- (44) 

Following the same procedure, in the Case R\ — Rlower mantle We obtain: 

^ « -0.096 (£--f). (45) 

From (|39| ) we find an upper bound on /: 

/ < 1.8 x 10" 4 . (46) 

Again, the condition of MAWE validity, Eq. (|38|) , is satisfied automatically with these /. 
Further, the upper limit on the mirror matter density p x becomes 

p x < 0.03 g/cm 3 . (47) 

Correspondingly, the bound on the mirror matter is 

M x 



< 3.8 x 10~ 3 (for R 1 = 0.89/?®). (48) 



Comparing our bounds ([EJ) and ( |4"7| ) with the limit (^) we see that the latter limit 
is significantly weaker than the former two (although the limit (|29|) has been obtained for 
the inner core, it is also valid for the outer core because of the similarity of their chemical 
composition; we also would not expect substantial changes of this limit in the case of the 



lower mantle). Therefore we conclude that Eq. (34), Eq. (f|4|), and Eq. ([|8|) represent our 
final upper bounds on the mirror matter mass. Being the largest of the three, Eq. (^) can 
also be considered as the most conservative, radius-independent upper bound on the mirror 
matter mass in the Earth. 



IV. DYNAMICAL MANIFESTATIONS OF MIRROR MATTER 

In this section we are going to analyze constraints on mirror matter that follow from the 
precision measurements of the Earth's gravitational field. Such constraints may arise if the 
mirror matter, for some reason, shifts away from the centre of the Earth. 

The motion of the mirror matter in the Earth would be controlled mainly by the Earth's 
gravity field; inside the core, in a first approximation this field grows linearly with the radius: 

g « kr, k = 3.6 x 10~ 6 s~ 2 . (49) 

(The net tidal force exerted on the mirror matter by the Moon and the Sun is negligibly 
small as will be discussed later.) Consequently, the period of the mirror matter motion 
inside the Earth T x would be independent of radius and given by 

T x = -JL = 3311 s w 55 min. (50) 
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This period gives us the time scale for the variation Sg of gravitational acceleration at the 
Earth's surface caused by the possible mirror matter motion. 

To determine the amplitude of the gravity variation suppose that the amplitude of the 
mirror matter motion is h. Then the amplitude of the gravity variation is given by 

The above equation holds exactly only for the circle on the Earth's surface which lies in the 
plane of the mirror matter motion. For points outside this circle Eq.(5T) can still be used for 



order-of-magnitude estimates. Requiring that the gravity variation should not exceed the 
observational limit, Sg/g ^5 1CT 9 (see e.g. jl2[0 we obtain an upper bound on the amplitude 
of mirror matter motion: 

, ^ „ /3.8 x l(T 3 M e \ 

h< 1.7xi — m. (52) 

What physical factors could lead to the off-centre shift of the mirror matter? Let us start 
by discussing the possible effect of Moon's gravity. The two key quantities to be considered 
are the tidal torque and the net tidal force exerted on the mirror matter by the Moon. 
The effect of the tidal torque would be to slow down the spinning of the mirror matter, in 
analogy with the ordinary tidal torque that brakes the Earth's rotation. The details of the 
effect depend on the poorly known characteristics of the mirror matter such as its angular 
velocity, elastic and dissipative properties etc.; we will not dwell on these. 

On the other hand, the net tidal force and the corresponding off-centre shift can be 
estimated without the knowledge of mirror matter properties. First of all we note that if 
the Earth was spherically symmetric then there would be no net tidal force acting on the 
mirror "ball" placed in the Earth's centre. However, if the oblateness of the Earth is taken 
into account then the net force at the centre is non-zero and thus the centre of the Earth is 
not an equilibrium position anymore. The new equilibrium position for the mirror matter 
can be found from the condition of balance between the net tidal force and the gravitational 
attraction of the mirror matter by the Earth. 

To find the new equilibrium position it is convenient first to find the point in the Earth 
where the lunar tidal force vanishes ("the tidal centre"). Due to oblateness of the Earth the 
positions of the Earth's centre of mass and the tidal centre are shifted relative to each other 
by a short distance b (hereafter we ignore the 18° inclination of the Moon's orbit relative to 
the Earth's equatorial plane): 

b = j 2 ^l ~ i7i m (53) 

where 

J 2 = ^—^ ~ 1.08 x 10- 3 (54) 

is the Earth's dynamical oblateness, R ~ 3.84 x 10 8 m is the distance between the Moon and 
the Earth, C and B are the moments of inertia of the Earth with respect to the principal 
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axes. Therefore the centre of mass of the mirror matter will move away from the Earth 
centre (and also away from the tidal centre) by the distance 

where m is the Moon's mass. Inserting Eq.(|53D into Eq.(|55D we obtain 



htt I ffi sa 8.1 x lO' 6 m. (56) 

The off-centre shift of mirror matter would create periodic variations of the gravity acceler- 
ation on the surface of Earth: 

*» s —Ml „ L3 x io-» * , (57) 

which is far beyond the observational limits. Thus we have shown that the effect of a net 
tidal force due to the Moon is negligible. A similar result holds for the solar effect: 

h « ° e » 10" 8 m 58 

km V ; 

where M « 2 x 10 30 kg is the solar mass, R ~ 1.5 x 10 11 m is the distance between the Sun 
and the Earth. As expected, the effect of the net tidal force due to the Sun is even smaller 
than the lunar effect. 

We now consider non-gravitational interactions that could possibly cause an off-centre 
shift of the mirror matter. Let us start by analysing the possible role of meteorites and 
meteor showers colliding with the Earth. Suppose that as a result of such a collision a 
momentum p is transferred to the Earth. Then the mirror matter (assumed to be at rest in 
the centre) would receive an initial velocity u = p/M @ relative to the Earth. Therefore, the 
off- centre displacement would be equal to 

h °" = 7Pr AfeTE' (59) 

What could be the magnitude of pi The maximal velocity of a Sun-bound colliding object, 
relative to the Earth, is v max ~ 73 km/s. The heaviest meteorite found on the Earth has 
the mass m max ~ 60 ton. Inserting these values into Eq. fl59|) we obtain: 

h col < mma * v ™ x « 3.8 x 10- 13 m. (60) 
M^Vk 

The variations of surface gravity acceleration caused by such displacements are many orders 
of magnitude beyond observational accuracy. Note that the impact of meteor showers would 
be much less than the estimate ( j60| ) since the total mass of even the most copious showers 
is significantly less than m max . 

In the case of still heavier meteorites which disperse after hitting the Earth the mass 
can be estimated only indirectly (see e.g. |13|). For instance, the meteorite that created the 



Arizona crater (with diameter of 1207 m and depth 174 m) had the estimated mass between 
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60 and 200 thousand tons. The mass of the Tunguska meteorite (fell in Siberia in 1908) 
was at least 1 million tons, its speed 30-40 km/s; if we insert these values into Eq. ( p9|) , we 
obtain 

h col ~ 10- 8 m; (61) 

corresponding values of Sg/g are still completely negligible (even in comparison with the 
variation of g due to the Moon's tidal effect, Eq. (|57|)). 

Let us now turn to another possible mechanism of the mirror matter motion. Can earth- 
quakes cause translational oscillations of mirror matter around the Earth centre? Large 
enough earthquakes are known to excite free vibrations of the Earth; the study of these vi- 
brations has become one of the most important pieces of information about the Earth interior 
(for more details see e.g. P,^4|-|T6|). These vibrations are classified into two categories: 

— toroidal, in which only shear strain is present so that density is not perturbed; they 
are denoted by r TJ; 

— spheroidal, where both shear and volume deformations arise, denoted by r 5j. The 
indexes r and I (for both S and T modes) have the same meaning as the radial and orbital 
quantum numbers of the hydrogen atom. 

Toroidal modes do not lead to gravity perturbations so they cannot excite oscillations of 
mirror matter. On the other hand, spheroidal modes might cause the excitation of mirror 
matter oscillations through the gravitational coupling. Note that we should distinguish be- 
tween two possible types of mirror matter oscillatons: a) bulk vibrations in which the centre 
of mass of mirror matter stays at rest in the Earth's centre and b) translational oscillations 
where the centre of mass of mirror matter oscillates around the Earth's centre in the Earth's 
gravitational field. We cannot say much about the spectrum of bulk vibrations without 
knowing the detailed structure of mirror matter (i.e., its density, elastic and dissipative 
properties etc.); for this reason we leave them out of our consideration. On the contrary, 
the period of translational oscillations can be found and is given by Eq.(j5C|). 

Our next task is to find out if there are any spheroidal Earth eigenmodes that could 
resonate with translational oscillations of mirror matter. Note that here we deal with the 
case of parametric resonance so we need to look for the eigenperiod Te = T\/2 « 1655 s 
rather than Te = T\ (the parametric resonance in the case Te = T\ is weaker than for 
T E = Tij 2). The closest such eigenmode is 0S4 with the period of || 

T( S A ) = 1545.6 s. (62) 

In the time-varying gravitational field of mode the frequency of translational mirror 
matter oscillations also becomes time dependent according to the law 



where 



uj 2 {t) = 0^(1 + a cos 7*), (63) 



2tt 2tt 5p 



Sp/p is the amplitude of density variation in the 0S4 mode. 
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The onset of parametric resonance is controlled by the quantity s called amplification 
index (see e.g. 0]): 



1 Ifauji 

2n— 



e = 7 — 2uji 



(65) 



If the amplification index is real then the oscillation amplitude grows with time as exp st. 
In the opposite case s 2 < parametric resonance does not occur. 
Using Eqs. ( |63| , |64| , |65D we can find s to be 



(6£ 



P 



(0.28) 



(66) 



which is clearly imaginary. Thus we conclude that the condition for a parametric resonance 
is not satisfied and consequently there is no amplification of translational oscillations of the 
mirror matter. 

It can be shown that the condition of parametric resonance with the o^-mode takes the 
following form: 



where 



2tt 



7 



3233.25 



is the frequency of the o^-mode. One can see that the condition 
there is no resonance with the oSVmode either. 



(67) 
(68) 

is not satisfied and 



V. CONCLUSION 

We have investigated in detail geophysical constraints on the possible admixture of mir- 
ror matter inside the Earth. To this purpose, a method has been developed based on the 
Preliminary Reference Earth Model — the "Standard Model" of the Earth which describes 
its internal structure derived from the geophysical data in a systematic and self-consistent 
manner. If the density of the mirror matter is given, our method allows one to compute 
changes in various quantities characterising the Earth (such as its mass, moment of inertia, 
frequencies of its normal modes etc.). Comparing the computed and observed values of these 
characteristics, we can obtain for the first time the direct upper bounds on the possible con- 
centration of the mirror matter in the Earth. In terms of the ratio of the mirror mass to the 
Earth mass these upper bounds range from 3.8 x 10~ 4 to 3.8 x 1CU 3 depending on the radius of 
the mirror matter ball. We then analyzed possible manifestations of mirror matter through 
the variations of the gravity acceleration on the Earth surface. These variations could arise 
as a result of an off-centre shift of the mirror matter due to several possible mechanisms 
such as lunar and solar tidal forces, meteorite impacts and earthquakes. Our estimates have 
shown that variations caused by these mechanisms are too small to be observed. 
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In this work we have been based on a standard premise that mirror matter interacts with 
ordinary matter only gravitationally []; we have not relied on any other specific assumptions 
about the mirror matter properties. Therefore our results are valid for other types of hy- 
pothetical matter coupled to ordinary matter by gravitation only; an example is shadow 
matter introduced in string theories. On the other hand, the use of equation of state and 
other macroscopic characteristic of mirror matter could lead to more severe constraints on 
the mirror mass inside the Earth. 
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Note that mirror matter can also couple to ordinary matter through photon — mirror-photon 
mixing . An analysis of this interesting possibility is beyond the scope of the present work. 
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